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Q^ ■ The presence of fluctuations and non-linear interactions can lead to scale dependence in the 

^0 ' parameters appearing in stochastic differential equations. Stochastic dynamics can be formulated 

I , in terms of functional integrals. In this paper we apply the heat kernel method to study the short 

*vj ' distance renormalizability of a stochastic (polynomial) reaction-diffusion equation with real additive 

noise. We calculate the one-loop effective action and its ultraviolet scale dependent divergences. We 
show that for white noise a polynomial reaction-diffusion equation is one-loop finite in d = and 
1-^ ' d = 1, and is one-loop renormalizable in d = 2 and d — Z space dimensions. We obtain the one-loop 

renormalization group equations and find they run with scale only in d = 2. 
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Many examples abound where particular spacetime distributions of matter are selected over a wide variety of 
seemingly possible choices. In many of these cases, these patterns are well described, and their temporal evolution 
'^ ', accurately modelled, by non-linear partial differential equations subject to noise, or equivalently, by stochastic partial 
P^ ■ differential equations (SPDEs). Concrete examples can be found in the domains of pattern formation, chemical chaos, 
O \ biological morphogenesis, and flame- front propagation, just to name a few [yj-^l. As argued in Ref. y the effective 
potential is a superb tool for studying the onset of pattern formation (i.e., symmetry breaking) about the static and 
spatially homogeneous solutions of the SPDE. This quantity takes into account both interactions and fluctuations to 
a given number of loops. (The loop expansion is a controlled expansion in the amplitude of the fluctuations.) To 
K*" ' go beyond the study of symmetry breaking and static homogeneous configurations, it is crucial to be able to include 
\l \ and account for time development and spatial inhomogeneities. Ultimately, we are interested in studying the late 
^1^ ■ time behavior of the solutions of SPDEs; we would like to know if there is an asymptotic steady state, or state of 
equilibrium, as this impinges on the late time behavior of the emergent pattern. The effective potential does not yield 
^~^ ' . this information, and one must turn to the effective action. The effective action contains all the dynamical solutions 
f^ ' of the SPDE and their asymptotic behavior. It obeys a variational principle wherein its first variation yields the 
^^ . SPDE. Its exact calculation, however, is very difficult. Along the way, one must first "regularize" its short distance 
^^ ' divergences that are dimension dependent. These must be identified, isolated, regulated, and if possible, renormalized 
2 ', by suitable redefinitions of the bare parameters appearing in the SPDE. The study of these dimension dependent 
divergences lends itself to a direct calculation of the renormalization group equations (RGEs) that encode the scale 
dependence (or "running") of these parameters. The RGE fixed points provide information about the asymptotic 
states. For these and other reasons, the divergent structure of the effective action warrants a detailed analysis in its 
Q ■ own right. 

^ , In this paper we study the stochastic reaction-diffusion (RD) equation by means of a functional integral represen- 

L| ' tation that will be used to calculate its one-loop effective action. In earlier work we calculated and analyzed the 

. ^H , effective potential for the RD equation (the effective action evaluated for special field configurations that are time 

^^ independent) Q. However, the effective potential does not provide information about the dynamics of the system or 

JH \ wavefunction renormalization. The effective potential can only signal the static homogeneous states around which one 

- - 1 can study the onset of pattern formation (i.e., symmetry breaking), but does not indicate which such state (if any) is 

dynamically accessible or most likely. To go beyond the limitations inherent in an effective potential analysis, and in 

order to judge the importance of the noise-induced symmetry broken ground state configuration, one should explore 

the space of dynamical (spatially inhomogencous) solutions of the system. For instance, are there solutions that 

indicate the system is "thermalizing" , or can we find steady state solutions? An important step in this direction is the 

analysis of the effective action undertaken in this paper. In this calculation we encounter short distance divergences 

that need to be regularized. We have chosen to carry out this regularization procedure by means of (generalized) 

heat kernel techniques. To do so, we will be required to calculate the first (integrated) Seeley-DeWitt coefficients. 
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In addition to the physical RD field itself, we must also calculate the contribution from the ghost field, a necessary 
ingredient in this formalism ^ . 

After the regularization step, we turn to the one-loop renormalization of the effective action. We find that for 
additive white noise RD systems are one-loop finite in d = and d — 1 space dimensions, and (at least) one-loop 
renormalizable in d = 2 and d — 3 space dimensions for polynomial reaction kinetics. There is no wavefunction 
renormalization in any of these dimensions, irrespective of the degree of the reaction polynomial Q. Moreover, the 
ultraviolet renormalizability requires including a bare tadpole Aq, or constant term, in the equation of motion. 

Application of heat kernel methods to stochastic field theories far from equilibrium may not be familiar, so a few 
comments regarding them may be useful. These techniques have been used primarily for calculating effective actions 
and the one- loop physics in quantum field theories (QFTs) and in curved space quantum field theory |0-Q. In the 
quantum domain, one is interested in computing quantities such as the effective action and the effective potential, 
which provide crucial information regarding the structure of the underlying theory at different length and time scales 
and are important in assessing the theory's renormalizability (or lack thereof), the determination of the running of 
couplings and parameters, patterns of spontaneous and dynamical symmetry breaking, and the structure of short 
distance (ultraviolet) and long distance (infrared) divergences [|l3|-|l7|. Moreover, for renormalizable theories, the 
computation of the effective action (actually, only its divergent part is needed) can be used to extract the RGBs that 
govern the scale dependence of the couplings and parameters appearing in the theory |1§H2C|] . Though perhaps better 
known in the context of these fields, these same techniques can be generalized and applied to reveal the corresponding 
one-loop physics associated with stochastic dynamic phenomena and to systems subject to fluctuations. 

As a key technical step, we need to obtain the Green functions for certain higher-order differential operators that 
are neither elliptic nor hyperbolic. (The operators treated here are second-order in physical time derivatives and 
fourth-order in space derivatives.) We set up a generalization of the Schwinger proper time asymptotic expansion 
to compute these Green functions and obtain explicit expressions for the first Seeley-DeWitt coefficients associated 
with these operators. We carry this out employing the minimal representation in which only physical and ghost field 
variables appear in the action [Q. 

At this stage we emphasize the following: We start from a given "classical" SPDE, together with a specification 
of the type of noise and its correlations. Our aim is to map this SPDE to an equivalent generating functional and 
effective action, at which point QFT techniques can be used. We emphasize that our approach and use of field theory 
is rather different from that initiated some years ago, whose aim is to derive the SPDE and the noise correlations 
from microphysics by mapping a classical master equation to "second-quantized" variables, and then finally to an 
action principle [pl|-p3[. 

We conclude with a summary and discussion of this work, and prospects for further development. 

II. STOCHASTIC FIELD THEORY FOR REACTION-DIFFUSION EQUATIONS 

In this paper we demonstrate that heat kernel techniques can be used to compute the ultraviolet divergent terms 
arising in the one-loop effective action associated with field theory formulations of SPDEs. In particular, we apply 
these methods to the class of single component reaction-diffusion equations 

D<l>{x)^V[cP{x)]+ii{x), (1) 

where 

= |-.V>, (2) 

and 

^ A- 

Here ri{x) is a noise term with normalized probability distribution given by P[r]]. We employ the shorthand notation 
X = {x,t). Furthermore, ¥[(/)] is a polynomial of degree N in the concentration (or field variable 0) and the Aj's 
are a set of reaction rates. For convenience we have placed the decay (or "mass" term) into the potential and have 
treated it as another coupling: Ai = —m?. It must be noted that this equation has the form of a purely dissipative 
system and has a bona-fide potential energy term V[4)] pq |. Thus, it makes sense to calculate an effective potential 
for constant field configurations, as well as an effective action for inhomogeneous fields. Both the effective action and 
the effective potential are derived by means of a field theory for this SPDE. 



Before continuing, we should point out that there are reasons to beheve that as a phenomenological equation, the 
RD being considered here, (and others structurally similar to it), might not be adequate to describe certain two-body 
annihilation processes, or pair reaction kinetics, since recent derivations based on master equations show that the 
SPDE in question should actually be complex with imaginary noise (leading to negative noise correlations) ||6|, |24| , p5| . 
On the other hand, for processes involving particle clustering, these same derivations yield real stochastic equations 
and noise, as well as positive noise correlations ||26| . Of course, there are many situations in which a microscopic 
derivation of the SPDE is entirely out of the question, cither because explicit knowledge of the microscopic details is 
lacking and/or because the random fluctuations owe to uncontrollable contingencies. In these situations the benefit 
of adopting a phenomenological strategy should be self-evident. Finally, the application of this equation need not be 
restricted to just chemical diffusion [p7| . 

For homogeneous and static concentrations it is sufficient to study the effective potential Q. In this paper we 
complement that analysis by making use of the effective action to consider inhomogeneous and time dependent field 
configurations. In the minimalist approach (see Ref. M) one starts with the normalized generating functional Z[J] 
encoding the stochastic dynamics described by the RD equation (|l|). This involves the RD scalar field (j) plus the 
Jacobian determinant (denoted here by J7) and its adjoint (JT"^), (these determinants arise from a change of variables) . 
The generating functional (partition function) is given by P,15[ 



Z[J] 



im exp{[-S[cj)] + /d"a; Jix)c^ix)]/A} VJj^ 



jm exp{-S[cb]/A} VJJ^ 
where the "classical" action is [valid for Gaussian noise: G,j{x, y) = {rj{x)r]{y)) = A 52 (a^, y)] 

Sm - I /d'^Ti fd"x2 {D(t>{x^) - F[^(a;i)]) g:^\x^,X2) (Dcj^ix^) - V[<f>{x2)]), 



(4) 



(5) 



with n = d + 1 and d the number of space dimensions (we will keep d as a free parameter throughout this paper). 
For a general SPDE there may be non- vanishing contributions from the "ghost" fields (Jacobian determinants). We 
follow here the discussion of Ref. Q to separate the noise two-point function into the product of a shape 32(2;, y) and 
a constant amplitude A. Irrespective of how we decide to normalize the shape, the constant amplitude A is always 
the loop- counting parameter of the perturbation expansion ^ . A loop-counting parameter is very useful in organizing 
such a perturbative expansion. Moreover, any symmetry that is present in the classical action (W is preserved at each 
order in the loop expansion since the loop-counting parameter multiplies the entire action (and the source term J) 
in (H). One of the advantages of the minimal representation is that it leads to this natural identification of the noise 
amplitude ||,|5|. 

We introduce the generating functional for connected correlation functions W[J] and its Legendre transform, the 
effective action r[(j}] |la] , (note the explicit factor of the noise amplitude A) 



with 



W[J]^+A logZ[J], 

r[0] = -w[j] + /d"x Jix) {4>[J]ix) ~ <^[o](x)} 



6W[J] 
SJ ' 



m = 



The barred fields (/)[J] and (f>[0] are the solutions of the equations of motion 






0[,/] 



Jix), 



and 






= 0, 



m 



(6a) 
(6b) 



(7) 



(8) 



respectively. It is usually assumed that the former equation has a unique solution 4>[J] (at least for small J), and that 
for vanishing source (J — 0) the unique solution is the vanishing mean field (j)[0] = 0, i.e., {(j){x, i))j=o = 0, where the 
angular brackets denote the stochastic average with respect to the noise probability distribution V [77] . This is valid 
for a symmetric ground state. We next expand the action (ra) about this background up to quadratic order in the 
stochastic fluctuation 54> = (p — cf). We can carry out a perturbation expansion in the small parameter A and compute 
the one-loop effective action to obtain [Of n 



^ We now drop the overbar on 4> with the understanding that this stands for the conjugate field of J and not the field appearing 
in the classical (zero-loop) action. 



r[0] = 5M - 5[0] + :^ (Trlog41[0] - Trlog^^'jjO] - \og J[4>] ~ logj^tf^] + iog^[o] + log jt[o]) + 0{A^) 

= 5[</.] - 5[0] + r(i-'-^)[0] + T%-J:f[cp] + 0(^2), (9) 

where the matrix elements of the Jacobi field operator 'S'^ddl^] ^^^ 

We have anticipated the appearance of divergences in the one-loop contribution to the effective action, arising from 
both the physical and ghost fields, and have supplied it with a cut-off e. 

The notation ^[O] is actually shorthand for S[cj)[J = 0]], and for a symmetric ground state one typically has 0[O] = 
and S[0] = 0, unless there is a "tadpole" in the classical action. In fact, when looking for mean field solutions 
of the zero-loop equation of motion, wc will find it convenient to consider a non-vanishing value of the mean field 

(2) 

(f)[0] = Wo 7^ and will study fiuctuations about this mean value. The terms involving S[0] and S'f;pj^[0] appear due 
to the normalization factor in (0). The notation "Tr" stands for the trace and indicates that we are to take the 
(time and space) coincidence limits X2 —^ x and xi — > x, followed by an integration over the common limit x. The 
one-loop effective action will contain divergent terms and it is precisely these terms we wish to isolate and compute. 
We have collected all such divergences into the expression F^ ~ °°^' and we regulate them by means of a cut-off e. The 
finite terms are collectively represented by ^finite''' ■ There may also be higher-loop contributions, denoted by 0{A'^), 
whenever we need to emphasize them explicitly. Although these latter contributions are important for constructing 
the full effective action, that calculation is beyond the scope of the present paper. 

In order to compute the one-loop effective action we need to obtain S^J^^ [0] . This Jacobi field operator is diagonal 
in coordinate space. For the purpose of this calculation and in the interest of simplicity, we consider the case of white 
noise. (Colored noise can be dealt with, but it brings in time and space derivatives of the shape function, which 
complicate the heat kernel analysis.) For white noise we have {'rj{x)'q{y)) — 2Do d"{x,y) and therefore we can write 

G^{x,y) = 2DnS"{x,y), =>A = 2Do, g2{x,y) = d"{x,y), and g^\x,y) = 6'' {x,y), (11) 

which fixes the noise normalization. 

The Jacobi field operator corresponding to the RD equation is easy to calculate starting from the classical action. 
We simplify notation and write the zero-noise action as 

S[cf\^\JA^x {D<f>-V[4>]f = ]^jd-x {dt(^-vV^<f>~V[^])\ (12) 

The Jacobi operator for the physical field, 5'£^j^((/>, xi, 3:2), is given by 

4eid('^:2;i,X2) = [{-dt - ^V2 - V'[<l,{x,)])(dt - ^V2 - V'[cp{x,)]) ~ V"[cp{x,)]{Dcp{x,) - V[<j,{x,)])\ S^ix^^x^), (13) 

where V'[<l^] = ™ and V"[cb] = ^. 

For the ghost field the Jacobi operator is given by |j] 

SSosti^^^l^^2) ^ {i-dt - UV' - V'l^ix^Mdt - UV^ - V^'[(/.(xi)])}<5"(.Ti,X2), (14) 

and its determinant can be written as 0] 

det[5gL(0,a;i,X2)] = J[4'] JH4>]- (15) 

In order to carry out the perturbation expansion we also need to consider the "free" case defined by the limit 
V[(t)] -> 

41(^1.^2) = [i-dt - lyV^) {dt - lyV^)] S^x,,x2) ^[~d^ + iiyVY] S'\x,,x2). (16) 

Free physical fields have the same Jacobi operator as free ghost fields, so that as V[(j>] — > 0, the physical and ghost 
field contributions to the effective action cancel. We now look ahead a little: as the Jacobi operator S^^^^{xi,X2) 



contains fourth order space derivatives (hi-harmonic operator)^ rather than second order derivatives, the behavior of 
the DeWitt-Schwinger expansion [[7|-[ll| is qualitatively different in that it includes fractional powers of the Schwinger 
proper time parameter. 

We now calculate the mean field wo (i.e., the background field) by studying the solutions of the classical equation 
of motion, which is given (for arbitrary source J) by 



5S\ 



> I 



i-dt - lyV^ - V'[(j)ix)]) {D,p{x) - V[(j){x)]) = J{x). 



0[j] 



If the source vanishes and the mean field is homogeneous and static, we have 

V'[v^] V[vo] - 0, 



(17) 



(18) 



which always has at least one real solution [g|. 

In order to calculate the one-loop effective action for RDs, one must include the contribution from the "ghost" 
fields. These "ghost" Jacobians are given by |B| 



5V 



and J^ ^Aet{D^ -^-^\ 



J ^Aet\D 
We can now complete the formal calculation of the one-loop effective action. We have 



T[4>\ - S[4>\ - S[vo] 



A 



Trlog4cid[ 



-Trlog(.t_^)(,_|, _(,..„) 



0{A% 



so the one-loop effective action receives one contribution from the physical field 



t(1 — ioop) 



A 



field o 

and a contribution from the ghost field 

-^{l — loop) -^ 



,(2 



(2) 



= IT Trlog5;,:;,[0] - Trlog5;,:i,[«o] 



r 



ghost 



Trlog(i?t_^U^_^) _Trlog(i?t 



-^Un-^ 



5(j) J 



5(f) 



(19) 



(20) 



(21) 



(22) 



We will soon see that individually, each contribution has complicated Seeley-DeWitt coefficients, but when taken 
together, the physical plus ghost sectors yield simple net Seeley-DeWitt coefficients. 



III. COMPUTING THE ONE-LOOP EFFECTIVE ACTION 

In this section we construct the regulated expression Fg ^ °°^' [cj)] for the RD equation. We follow closely the 
DeWitt-Schwinger (DS) proper time formalism to analyze the ultraviolet divergences |7[- |ll| , [l4| -[l6|| . (We have striven 
to keep this section self-contained.) 

In this formalism the integral representation for T^^^^°°P'[(j}], eq. (Q), involves a fictitious "time" parameter s 
(denoted as Schwinger proper time). To this end, we define the following function, where A is any operator 



go,{A) = / ds s'^-^e-"^ = A"" r(a), 
Jq 



with r(Q;) the Gamma function. We consider the limit a — > 

1 



ga{A)^ --7-logA, 
a 



(23) 



(24) 



where 7 = 0.577..., is Euler's constant. Although this integral is divergent, the difference of two such integrals is finite 
and well defined 



/■+°° ds 
\im[g^{B) - g^{A)] - log A - logS = - / — ( 



,-sA_^-sB-^^ 



(25) 



and comparing with (^, the desired proper time integral for the one-loop effective action is given by 



-p{l — loop) 



[</-] 



-^{l — loop) 



eliost 



[<!>] 



A [+°° As f ,„ ,,„. \ A [+°^ ds 



where the "Hamiltonians" in the exponentials are the Jacobi operators 

:.(2) 



j "!^Tr(e-^^«"--e-'^[^°U-), 



Afield - S\^J,MJ]] - {D' - V'){D - V) - V"{D<P ~ V\ 



(2) 



[^oifieid - stiAm\ 



H,i,ost = S'Zjcj^iJ]] = (Z?T _ v'){D - V), 



[Ho 



J ghost 



s: 



(2) 



ghost L'/'[0]], 



(26) 
(27) 



(28a) 
(28b) 
(28c) 
(28d) 
(28e) 



as can be seen by comparing (Bq) with (m, ([l3|), (M), and (Iq). To proceed with the calculation, we need an explicit 
form for the operators e~'^^ (that is, for 6^"^^'=''', e^^I^^lfioid, g^'^^si'™', and e^*' olghoat)^ or rather, their matrix 
elements, so that we can take the indicated traces. To solve for them, we note that e^'^^ is the exact solution of the 
following operator differential equation 



OS 



(29) 



If one takes matrix elements in the spacetime coordinate basis |a;) = \x, t), inserts a complete set of states, and makes 

(2) 

use of the diagonality of H in the coordinate basis [note that S'fjpjjj [(/>] is proportional to (5"(xi,a;2)], we obtain 



Hix) (x|e-^^|x') == / dy{x\H\y){y\e-^"\x') = {x\H e-'"\x') = -yjx\e-'"\x') 



(30) 



or equivalently 



H(x) G(x,x'\s) ^- — G(x,x'\s), with G(x,x'\s) = (xle 

OS 



-'"\x') 



(31) 



This latter equation defines the Green function G{x, x'\s) in terms of the matrix element of the operator e^'* in the 
coordinate representation. These steps can be repeated for the other Hamiltonians. Fortunately, for the purposes of 
the present work, it is not necessary to solve this equation exactly (for either -fffidd or -ffghost)i as we are interested in 
the short distance divergent part of the one- loop effective action. What we will do instead is solve the "heat" equations 
(p9|), (pO|), and (pi]) adiabatically by expanding in small positive fractional powers of the proper time variable s, which 
is where all the ultraviolet divergences are to be found (different techniques are required if one is interested in the 
infrared limit). Nevertheless, to get "off the ground" it will be most useful to have the exact solution to ( POJ ) in the 
free limit (V^[(/)] ^^0). We now turn to this task, which entails solving exactly (p9) with /if free- 

Since equation ( pO|) looks like a heat equation in a n + 1 dimensional spacetime (parabolic equation in the proper 
time variable), we know how to solve it together with specified boundary and/or initial conditions. In the free field 
limit (G — > Gfree), we must solve the following equation 



-df + {iyV 



2\2 



_9 



Gf,,,{x,x'\s)^6''{x-x')6{s), 



(32) 



subject to the boundary (initial) condition Gfreo(a;i2;'|0) = S"{x — x'). Strictly speaking, the Green function depends 
on both arguments x and x' , but due to the translational invariance of the dynamical equation, we have G{x, x'\s) — 
G{x — x'\s) = G{x — x',t — t'\s) and it suffices to treat it as a function of one spacetime coordinate. We can always 
restore its dependence on both spacetime arguments at any time. 

The formal solution for Gfroc, expanded to fourth order in x — a;' (along the diagonal), is given by 



Gfrco(a;,a;'|s) = Ad exp 



4s 



1- 



C± 



Cd.2 (|x-ff)2 



2d 



M{d + 2) 



O 



V'- s 



■l\2\Z 



,,3/2 



, (33) 



where 

U»/ 2(21)'' r(f) l,rV * T(i) ^ ' 

Details of the calculation leading to the final expression for G^^''^"'' {x , x'\s) are given in Appendix ^. 

We also wish to point out that for static and homogeneous background fields (vo) the computation of Gfioidi^o] ^^nd 
Gghost[i'o] is not much more complicated than that for Gfrcc- Details are presented in Appendix |B[ These static and 
homogeneous calculations allow one to compare with the effective potential formalism of B , and serve as a check on 
the current effective action calculation. 

We adopt the following ansatz to perturbatively solve the heat equations (|30|) for small s (adiabatic approximation) 

Gficid(a:, x'\s) = Gfroo(a:, x'\s) /ficid(a:, x'\s), (35a) 

Gghost(a:,a:'|s) = G{rccix,x'\s) fghostix,x'\s), (35b) 

where 

+ 00 

fficid{x,x'\s) = ^s5 b^{x,x') = &0 + s^bi + sbi+ 0{s^), (36a) 

+ C50 

/ghost(a;, x'\s) = 2_] ■s^ ai{x, x') = aq + s^ai + s ai + 0{s^), (36b) 

1=0 

are asymptotic series in half- integer powers of the proper time with coefficient functions bi_ and ai (called "Seeley- 
DeWitt" coefficients) . Note that we have had to consider fractional powers in this small s expansion. [By considering 
simple cases it is easy to convince oneself that for a differential operator of order n the "heat kernel" expansion should 
start with an overall factor proportional to s~ '" and then contain subdominant terms that are integer powers of 

In principle, these coefficients can be calculated to arbitrarily high order by solving recursion relations obtained 



by substituting (35a) and (|35q) into (pOf) for -ff field and if ghost, respectively. The boundary (initial) conditions 



Gfie\dix,x'\0) = Gghost(2:i2:'|0) — (5"(x7P) imply that 5o = f and Oq = f. These coefficients start the Seeley 



DeWitt hierarchy and allow for a complete determination of the Seeley-DeWitt coefficients appearing in ( |36a| )-(36b). 
For second-order differential operators this procedure has now become automated ^8| . For fourth-order differential 
operators considerably less is known |29]. 

In practice, we will see that only the first Seeley-DeWitt coefficients are germane to the problem and that it is 
sufficient to find the "integrated" Seeley-DeWitt coefficients. This permits us to calculate the relevant coefficients by 
means of a technique based on the Feynman-Hellman formula |30[ , which can itself be viewed as a specialization of 
the Baker-Campbell-Hausdorff formula |pl| . 

We regulate the one-loop effective action by cutting off the lower limit of the proper time integral 

p(i-ioop)[^] = -^ T" ^ Tr (e-'^^-- - e-^I^"!--) + ^ T^ ^ Tr (g-^^^-- - e-^[^°l«--) , (37) 

where we can identify e = f/ri^^^^^jj and Te °°^'[(/)] with F^^ °°^ [0]. As the product sH must be dimensionless, we 
deduce that s has engineering dimensions of (time)^ or equivalently, (frequency)^^. In this stochastic field theory the 
cut-off can be taken to be a frequency scale i^cut-off , and this identification allows us to compare between these two 
types of cut-off (proper time versus frequency). Since these theories are not Lorentz invariant, a frequency cut-off is 
not "quite" interchangeable with a wavenumber cut-off (more on this point below) . 

Substituting the above ansatz (35a)-(35b) into (pTh, making use of the explicit form for Gftoc (p^), and expanding 



out the first terms of /field and /ghost yields 

p(l-ioop)j^j _ 1 f+°°ds 



A 2 J, s 

1 r+°^ ds 

2 



2^' 



d"x [G{rcc{x,x'\s)] ([/fieid(2;,a:'|s)] - [/ghost (a;, x' I s)]) 

d"x [Gfree(a;,a;'|s)] 
s2 [bi{x,x')] - si [ai{x,x')] + s [bi{x,x')] - s [ai{x,x')] + 0{si)] 

"■x (si [ci{x,x')] + s[ci{x,x')]+0{si)y (38) 



ds 



In arriving at this last expression we have used the fact that the coincidence hmit of the free heat kernel is 

[Gt,,,ix,x'\s)]^AdS-i-^, (39) 



which follows immediately from (Bq). We have also made use of the standard notation to express coincidence limits. 
Given any function h(x,x'), we write 

lim h{x,x') = [h{x,x')]. (40) 

x' — ^x 

(Although we also employ the square brackets to denote arguments of functionals and functions, and to group expres- 
sions, the intended meaning should be clear from context and there should be no confusion.) We have defined net 
Seeley-DeWitt coefficients 

ci =bi ~ai, yi = 0,1,2,.... (41) 

2 2 2 ^ ' 

The Seeley-DeWitt coefficients bi, ai, and ci are functions of the mean field (j){x,t) and its derivatives, and as 
remarked above, can in principle be determined by solving a recursion relation resulting from inserting the ansatz 



(35a) and (35b) into (30). However, to obtain the form of the divergences of the one- loop effective action we need not 
evaluate these coefficients. It suffices to calculate the lower bound (s -^ 0) of the proper time integral. In the limit 
e ^ we find that the divergent terms in the RD effective action are given by 

r(i-'-rt[^] ^ ^^A,A n e-i Jd-x [c.{x,x')] -^Z^Jd-x [c^{x,x')] 




^^ „. , d"x[c,{x,x')]+ •••). (42) 

We now list the divergences in the RD one-loop effective action for the following space dimensions 

d = r(i-'°°P)[(/)] = -2Ao A log(f^2 ^) f ^^ j^^i^ (43^) 

d=l ri'^-^''°P^[(t>]^-2AiAe-i f dxdt[ci], (43b) 

d=2 r(i-'°°P)[0] = -^ A (2e-i I d^xdt [ci] -log(172 e) j d^x dt [ci]\ , (43c) 

d = 3 r(i-'°°P) [(/)] = -^ ^ ('^e-i /"d3fdi[ci]+4e-3 /"d3fdt[ci]y (43d) 

In all of these cases we only need to solve for the first two adiabatic (Seeley-DeWitt) coefficients ci , and ci; indeed, 
it is only the spacetime integrated net coefficients that are needed. In higher space dimensions, additional ci_ 's would 
be required. However, for most practical applications it is enough to consider Q < d < Z. (Moreover, earlier work 
regarding the one- loop effective potential for RD indicates that this field theory is non-renormalizable for d > 4 [g[). 
This dimension range covers the spatially homogeneous limit {d = 0), one-dimensional (linear) systems {d = 1), 
surfaces (d = 2), and bulk systems (volumes) (d = 3). We see that the divergences are of two types: (fractional) 
powers of the cut-off and logarithms of the cut-off. Only the latter can yield one-loop renormalization group beta 
functions and associated RGEs. 

IV. CALCULATION OF THE SEELEY-DEWITT COEFFICIENTS 

In this section we present a formalism that can in principle yield all the Seeley-DeWitt coefficients. As we have seen 
in the previous section, the calculation of the one-loop effective action involves solving the auxiliary partial differential 
equations of parabolic type (denoted as heat equations, even if the diffusion is non-standard). 

In this formalism (see Appendix O) the first step is to write 

Tr[exp(-siJ)] = Tr{cxp[-s(iJft.ec + 6H)]}, (44) 



where SH is a lower-order differential operator when compared to H or -ffftcc- We now apply a version of the 
Feynman-Hellman formula [ pOUSlI , as discussed in Appendix |y, to obtain 

Tr[exp(A + eB)] = Tr[cxp(A)] + eTr[S exp(A)] + — d^ Tr{B exp(M) B exp[(l - £)A])} + O(e^). (45) 

2 Jo 

This equation will be the basis for extracting the first two integrated Seeley-DeWitt coefficients. 
The second step is to realise that we only need to look at the difference 

Tr[exp(-s7Jfieid)] - Tr[exp(-siJghost)], (46) 

and write 

Tr[exp(-s7?gho.t)] = Tr{exp[-s(iffrcc + SHi)]}, (47a) 

lV[exp(-sHfieid)] = Tr{exp[-s(i/frec + SH^ + 6H2)]}. (47b) 

If we take the difference of the previous operators, the 0(e") term automatically cancels, as does the 0[e^{SHi)^] 
term, leaving 

Tr[exp(-siJficid)] - Tr[exp(-si/ghost)] = -esTr[(5i?2 exp(-si/freo)] 

+ -s^ / de Ti {SH2 exp(-£sFfrcc) SH2 exp[(l - ^)siJf,ec])} 
^ Jo 

+eh^ f di Tr {SHi exp{-isHf,,,) 6H2 exp[(l - £)si/frec])} 
Jo 

+0{e^). (48) 

We now make use of the explicit form of these "Hamiltonians" . We write 

TJghost = {D^ - V'){D - V) = D^D - (ZJt - V')V' + 2iyWV' ■ V, (49) 

where in the last term of the right hand side both V's act on everything to the right. We know that the free 
Hamiltonian is given by 

Hf,,, = D^D = [-d^ + iy^(yY], (50) 

so that we can identify 6H1 as 

SHi^ ^[{D^ -V')V' -2iyVV' -V]. (51) 

(Note that 6 Hi is a linear differential operator.) From the definition of the ghost Hamiltonian 

infield = i^ghost - V'iD^ - V), (52) 

we deduce the following form for SH2 

5H2 = -V"{D(I)~V). (53) 

(Note that 5H2 is a function, not a differential operator.) Now consider the first order perturbation [the 0(e) term] 

Xi = s 1:t{[V"{D4> - V)] cxp(-si/frec)}. (54) 



From the known form of the free kernel, [see, e.g., equation ( p3| ) or equation (All)], and the fact that 5H2 is a 
function, this reduces to 



1 d 



Xi = Ad s-5- 4 / d"a; \s V"{D(j) - V)] . (55) 



This implies that the first-order perturbation does not contribute to the Seeley-DeWitt coefficient c 1 , though it does 
contribute to ci. In fact, we can write 



d"x[ci] = / A''xV"{D(l)-V) + 



(56) 

This is actually the only contribution to the relevant Seeley-DeWitt coefficients. (There might have been additional 
contributions coming from those portions of the second-order term X2 that have space derivatives; fortunately they 
vanish, as we now verify.) Let us consider 

X2 = y / d£ Tr { [V"{D<I) ~ V)] exp(-s I iJ^ee) [V"{Dcj) - V)] exp[-s (1 - £) Hf,,,] } (57) 

+ s^ f d£ Tr { [V'iD^ - V)] exp(-s £ Ht,,,) [{D^ - V')V' - 2iy{W^ + (yV) ■ V)] cxp[-s (1 - £) H{,,,] }. 
Jo 

We can have any of the following cases 

- No gradients hit the free kernel: the term containing two factors of [V"{D(I) — V)] is proportional to s^ and can 
only contribute to C2, which is not needed in the present context. 



- One gradient hits one kernel: from equation ( |33| ) or equation ( All ), one can see that there is a factor of 
(x — x')i/{iyy/s) of order s'^". Such a term is odd under the interchange of x and x' and will vanish when taking the 
spacetime trace. 

- Two gradients hit the same free kernel: there will be contributions of the type 



Cd.i — 1= + {Cd 



cPs 



(58) 



which, after tracing with the free kernel, yield contributions proportional to s"^'^. Therefore, these terms contribute 

to C3 , which is not needed. 
^ . . . . 
Continuing in this manner, it is easy to convince oneself that there are no additional contributions to the required 

coefficients ci and ci. We can finally write 



/"d"a;[ci] =0, 

[ d"x [ci] = / d"x [V"{D(t) - V)]. 



(59a) 
(59b) 



Note that the present calculation only yields the integrated on-diagonal (x = x') Seeley-DeWitt coefficients and is 
insensitive to any term that vanishes upon integration. 

With a little more work along these lines, it is also possible to obtain the Seeley-DeWitt coefficients: [ai] and [ai] 
for the ghost field, and [61] and [&i] for the physical field. We only quote the results here 



I <rx[ai_\ = f d"x2Cd.i V, 

/"d"x[ai] = /"d"a; 

f d"x[bi] = f d''x2Cd,iV', 



^ {V'f + {D^ - V')V' 



d"x [61] 



d"x 



{V'f + (i?t _ y')y' + v"(r)<\> - V) 



(60a) 
(60b) 
(60c) 
(60d) 



V. ONE-LOOP RENORMALIZATION 



We have already calculated the (regularized) one-loop effective action ( |20D for the RD equation, and thus, we 
may now explore the renormalizability of this field theory following the prescription reviewed in pG] . In order to do 

so we must analyze the divergences of the one-loop effective action V\ ° . We must also keep in mind the fact 
that the bare theory (i.e., defined by the action (O)) does not depend on the arbitrary scale /i introduced by the 
renormalization scheme. Therefore, just as for the case of QFTs [18 19 1, we will derive a set of equations that govern 
the scale dependence of the parameters appearing in the RD effective action from the identity 
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d/i d/^ d/x 

where the O(yl^) indicates there wiU be higher- loop contributions to the effective action. In quantum field theory this 
identity does yield the one-loop RGEs since equation (pfl) can be expressed in terms of a sum of independent field 
operators (operator basis) and each coefficient of an element of this basis determines an independent RGE. 

As we have already calculated the relevant Seeley-DeWitt coefficients for the RD equation, we now turn to in- 
vestigate the one-loop renormalizability of its stochastic field theory. The renormalizability criteria are based on 
the following definitions. For renormalizable and super-renormalizable theories the counterterms needed to cancel 
the divergences are equal to, or fewer in number than the terms appearing in the zero-loop action, which implies 
that the Seeley-DeWitt coefficients are expandable in terms of the same operator basis appearing in the classical 
action. In particular, this basis set consists of {dt(t>, V^(^, l,(f>,(jP,- ■ ■ , 0^}. For non-renormalizable theories this cri- 
terion fails. That is, there are terms in the integrated Seeley-DeWitt coefficients that do not appear in the classical 
action ||jl5[|3|,P 



By comparing the zero-loop action (12) with the divergent terms of the one-loop effective action, it is easy to see 
that the latter do not contain any field operators not already present in the bare (classical) action. The divergent 
contributions to the one-loop effective action in d-dimensions are given by 

r(i-'°°^)[0; vo] ^^Ad -^^ fd^x di [ci(0) - c^M] + 0{e^-i) 

^Y^'^tSt /d'a;di{^"(0)[^0-n0)]+^"(«o)n«o)} + O(e'"^). (62) 

Some remarks are in order. First of all we point out the fact that the one-half Seeley-DeWitt coefficients of the field 
and ghost mutually cancel out. This cancellation is special to the RD system and does not take place for generic 
SPDEs. Secondly, the ill-defined quantity "e°/0" arising in d = 2 must be replaced by log(ri^e) — log(rj^/r2^^j_Qg). 
The dimensionfuU (but arbitrary) parameter Q is required to make the argument of the logarithm dimensionless. It is 
often more convenient to introduce a cut-off in wavenumber, rather than in frequency. In Lorentz invariant theories 
(QFTs, for example) these are essentially equivalent and it is usual to adopt units where the speed of light is one. 
In the RD system this would be inappropriate, since the equation is not Lorentz invariant. Instead, one notes that 
dimensionally 

[e] — [proper time s] — [physical time t]'^ = [t/]~^ [distance]**, (63) 

and therefore, in terms of a wavenumber cut-off A and a wavenumber subtraction point /i 

"670" ^ log(fi2^) = log(f}V^™t-off) = log(/iVA^). (64) 

(This observation is important when comparing different regularization schemes; for instance the effective potential 
calculation of ||] uses a wavenumber cut-off.) 

Thirdly, it is of great importance to study the type of divergence arising in the one-loop effective action for the RD 
equation, i.e., logarithms versus (fractional) power. From the above we see that the type of divergence depends on 
the number of space dimensions d. If d is odd, there will never be logarithmic divergences to one-loop order in the RD 
field theory; to get a logarithm, the space dimensionality must be even. A similar feature holds also for the one-loop 
effective action for QFTs in odd spacetime dimensions [p3,p3[ . We can conclude that the appearance of logarithmic 
divergences for specific space dimensions is not an artifact of the RD field theory, nor of SPDEs, nor is it an artifact of 
the regularization scheme we have employed. In QFT the RGEs yield the running of the coupling constants, i.e., give 
the scale dependence, if and only if, there are logarithmic divergences in the effective action. Thus, we can already 
predict that the parameters in the RD equation will not run in the ultraviolet region (to one-loop order) for odd space 
dimensions [ p3[ . 

Nevertheless, the one-loop effective action in odd space dimensions still contains divergences (though not for d = 1), 
which must be subtracted by suitable counterterms, but once this subtraction has been performed, the remaining finite 
part Y finite^ is independent of the subtraction scale /u. 

VI. RENORMALIZATION OF THE RD EFFECTIVE ACTION 

In this section we calculate the counterterms needed to renormalize the one-loop effective action. We first start 
with the bare classical action (|lj) for the reaction-diffusion equation 
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S[cf>] 



d"x{dt(t)-i'W^(t)-V[(l)]y 



d"2; df 



vV' 




The bare parameters (no subscript) are related to the renormahzed ones (denoted by a subscript R) as foUows 



A, 



z-^'\ 



\f + 5\j, 



with Z =l + 5Z, 



^ = vr + Sv, 



(65) 



(66) 
(67) 
(68) 



with i5Z, SXj, and 6i> the corresponding counterternis for Z, Aj, and v, respectively. (Our convention for the definition 
of the wavefunction renormalization constant Z does not follow the standard one in QFT |14,0|.) 

Our task consists in demonstrating that all the divergences appearing in the regulated one-loop effective action 
can be cancelled by suitable choices for these counterterms. In effect, we absorb the divergences into the (bare) 
parameters of the RD equation by renormalizing these parameters. If we write the bare action in terms of the 
renormahzed parameters and counterterms, and keep up to linear order in the counterterms, (which is sufficient for 
a one- loop analysis), we find 



S[ 



1 /■ / ^ A^ \ A7 /■ / ^ A-" ^ 



d"x dt(l)R - UrV' 



•}R 




N 



5v\I' 



OR 






.]'■ 



which can be written in a more compact and transparent notation as follows 



(5A,- 



Af(l-j) 
-5Z^— — 



S\<h\ = 



^jd"x{DRC^R - VR[<j>R]f -^jd"x{DRCl>R - VrIcJ^r])' 

I d"a; {Dr<j,r - Vr[M) 5v V^^r + f^ " 



SXi + 5Z 



Af(l-J) 



= SR[4'R,fi] + Ss[(f>R], 

where we have introduced the scale dependent renormahzed action SR[(f)R, /i] 



S 



R[<PR 



,^A 



d"x Z5 {fx)dtcj> - M^)Z- ip)V^<f> - V[Z-2cp, Af (m)] 



with 



N XR(n) 



3=0 



J'- 



(69) 
(70) 

(71) 
(72) 



the renormahzed (scale dependent) potential. 

The meaning of Dr is clear from inspection. The final equality in ^u) defines the finite, renormahzed action Sr 
and the divergent (but regulated) counterterm action Sg- The individual counterterms appearing in Sg will be used 
to cancel off the divergences arising in (62). We carry out this cancellation separately in each space dimension since 



each case leads to structurally different divergences [see equation (|6: 



A. d = counterterms and renormalization 



The case d = is very simple: in zero space dimensions there is no diffusion. There is a brief discussion in 
reference jS) and we do not belabor the point here except to reiterate that in d = the RD equation is one-loop 
renormalizable and one-loop finite. 
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B. d = 1 counterterms and renormalization 

In one space dimension the (formally) divergent effective action is given by 

r^'"'°°^^M = -2^Aie-i/^ fdxdt [ci]+0{A''), (73) 

with Ai = r(l/4)/[87r(7ri/j^)^/^]. We did not explicitly write this term in equation (|6|) since it vanishes identically. 

From our previous calculation of the Seeley-DeWitt coefficients we know that [ci] — (in all dimensions), which 
tells us that in one space dimension there are no divergences, that is, the theory is one-loop finite and there is no need 
to introduce counterterms. Since no renormalization is required, there will be no scale dependence in the parameters 
appearing in the RD equation. The beta functions, f3o = udO/dii, (encoding the scale dependence of the parameters) 
are therefore zero [at least up to order 0(^^)], and we have 

Z=1 + 0{A'^), (74) 

iy = i^R + 0{A'^), (75) 

X,^Xf + 0{A'). (76) 



C. d = 2 counterterms and renormalization 

In two space dimensions the divergent effective action is given by equation (|6S 

^AM - +^log(///A*)y'd2f dt [ci], (77) 

where A2 — l/(167ri'ij). 

From the calculation of the Seeley-DeWitt coefficients we know that for d = 2 we have 

[ci] = 0, (78a) 

[ci] = Vli[4>R]{Dq^R - VR[ct>R]), (78b) 

where we have written the Seeley-DeWitt coefficient [ci] in terms of the renormalized parameters as we are only 
working to one-loop order. Therefore, for the divergent part of effective action we can write 

r(i-'-rf [M = ^ log(M/A) J d^x dt VII[M{Db.^b^ - VrIM]). (79) 

In order to determine the value of the counterterms and to cancel them off, we must set 

r(i-'°°^)[0;v.] = -Ss[M + finite. (80) 

This cancellation can be made up to a residual finite but scale dependent logarithm. This is because the difference 
of two divergent logarithms can be finite and non-zero. The counterterms are proportional to log(/i/A), where /i is 
an arbitrary scale needed to render the argument dimensionless, but this scale need not coincide with fiQ, the other 
arbitrary scale needed to render the argument of the other logarithm, log(/io/A)i dimensionless 133]. 

If we perform this cancellation, we obtain the following //-dependent family of solutions for the counterterms SXj 

^A, log(MVA^) vm = ^A, log(^VA^) E ^^(^ ~ 1) ^''" -E-^^- (81) 

j=0 ^' J=0 ^' 

As we are working to one-loop order, we can set Z(fj,) equal to one in V^'. We can then read off the individual 
counterterms from this equation, using the linear independence of the basis {dtcf), V^(/>, 1, cj), cf)^ , • • • , 4> }, to obtain 
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SZ = + O{A'^), (82a) 

Siy = + O{A'^), (82b) 

SXo = ^ log(M/A)Af + 0(^2)^ (82c) 

-^Ai = ^ log(/x/A)Af + 0{A% (82d) 

<5A^_2 = ^ log(M/A)A;« + 0(^2^^ (82e) 

SXn-i=0 + O{A^), (82f) 

(5AAr = + O(yl2). (82g) 

In particular, we see that there is no wavefunction renormahzation nor diffusion constant renormahzation in two 
dimensions at one-loop order. The couplings associated with the highest and next-to-highest powers of the field 
(A^_i, A^) do not require renormahzation to this order. 

As pointed out above, due to the logarithmic divergence in two dimensions, when we subtract the divergences from 
the counterterm action, we are left with a finite /i-dependent logarithmic piece in addition to the renormalized action, 
that is 

r[0] = 5[0] + r(i-'-p)[</)] + r(^Jr)[0] 

= Sr[cI>r, A^] + ^ log (^^^ I d^S dt {Dr4>r - VR[<j,R])Vi^[(pR] + 0{A^). (83) 



We insert this expression into (61) to obtain the one- loop RGE 



In arriving at this equation, we have cancelled an overall common factor of the classical equation of motion, since 
Dr(J)r - Vr[(J)r\ 7^ in general. 

By collecting up the coefficients of the linearly independent terms in (|4|) we find the corresponding one-loop RGEs 
and beta functions in d = 2 to be given by 

|3z^^i^ = Q + 0{A^), (85a) 

/?, = ^^ = + 0(^2)^ (85b) 

/5ao = M^ = -F^A? + 0{A% (85c) 



d/x ^t^vr 

dAf _ A 
djj, SirvR 



(3,^ = ^^ = --^Xi + 0{A'), (85d) 



d/i 87rz/fl 



/3a„_. = M— f^ = -z ^N + OiA^), (85e) 



dA^ 
/3a„_,=M^^ = + 0(^2^, (85f) 

/3a„=M^ = + 0(^2). (85g) 

d/i 

Since there is no wavefunction nor diffusion constant renormahzation the set of one-loop RGEs can be summarized 
in the following way 
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A^^^ = -^^'[0«]+O(^^). (86) 

This equation agrees with the computation of the RGEs based on the effective potential, which was calculated in 
Ref. dl. [See equation (51) of that paper.] Furthermore, if we define fj, = /zoexp(r), the previous RGB becomes 



dyfl[0fl] A d^VR[(j>R] 



dr Sttvh d(/) 



+ OiA^), (87) 



which implies the fact that the one-loop RGB in d = 2 behaves like an anti-diffusion process in field space. 

At this point it is interesting to compare our results with independent renormalization group results obtained, for 
example, by Cardy in Ref. g. If a path integral is derived (along the lines given in ||2l|-p^) for the two-body process 
A + A —* inert, then the corresponding RD equation turns out to be given by Q 

D(j)=-2X(j)^ + ri{x), (88) 

where however, the noise must be pure imaginary. A renormalization group analysis of equation (Bq) shows that the 
field (j) does not require wavefunction renormalization, nor does the diffusion constant i' renormalizc. Our one-loop 
heat kernel computation performed for an arbitrary reaction polynomial, (85a, 85b), is in complete accord with these 
results (even though we treat real noise). Returning to (p8|), the only non- vanishing renormalization is that of the 
coupling A. It turns out that the one-loop renormalization group beta function is exact, and when expressed in terms 
of the dimensionless renormalized coupling gji is given by p] 

/3(5fl) = bgl, (89) 

in d — 2 dimensions, where 6 is a positive constant (the value of this constant is not specified in H). 
In order to compare these results, we define the following dimensionless couplings 

5.-^^, 0<j<N-2, (90) 

provided, of course, that for a given j the coupling constant Xf ^ 0. This definition together with the hierarchy of 
beta functions given in (85) show that the dimensionless coupling constants gj satisfy the following one-loop RGB 

where the overdot is shorthand notation for /id/d/i. We now consider an RD equation of the type given in ([l|) with 
real noise and for a degree-two {N = 2) reaction polynomial (0) V[<l)] = Aq -I- ^0^- This particular choice is made in 



order to be able to treat an RD equation as close as possible in structure to the one in (88). Apart from the imaginary 
versus real noise, the essential difference lies in the fact that we (must) have a tadpole term, whereas (IsS) lacks such 
a term. In this case there is only one dimensionless coupling which can be defined, namely 

A X^ 

and equation ( |9l| ) implies the following one-loop beta function for this dimensionless coupling 

f3(9o) = 9o{^- ^) = -9o{^) = 9l (93) 

This follows from (85c) and from the fact that Af ex A|^ = 0. 

Thus, for the purposes of renormalization and calculating the RGEs, this example demonstrates that it is equivalent 
to start from a complex or real SPDE, and that the field theory can be derived from microscopic principles or obtained 
by means of the procedure outlined in B . 
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D. d = 3 counterterms and renormalization 

In three space dimensions the divergent effective action is given by equation (p3) 



r(i-'°°p)[0^] = -2AA^t~^l^ f d^xdt [ci], (94) 

with ^3 = r(3/4)/[167r(7riyfl)3/2]. 

The net Seeley-DeWitt coefficient is given by 

[ci]^V"[MiD^R-VR[M)- (95) 

The vanishing of the index one-half Seeley-DeWitt coefficient means that the divergent effective action in d = 3 
becomes 

ri'-'^^P^M = -2AA^e-^l^ [d^xdt V'iMiDncl^n - V[<Pr]) + 0{A^). (96) 



In order to determine the value of the counterterms we must once again set 

r(i-'°°P' [M = -^,[0fl]+finite. (97) 

The last identification yields the following (/i-independent) set of counterterms 

6Z = + OiA^), (98a) 

Si^^Q + 0{A^), (98b) 

6Xa = -^AAsX^e-^/^ + OiA''), (98c) 

^Ai = -AAsX^e-^/^ + OiA''), (98d) 

SXn-2 = -AA^X^e-^^^ + OiA''), (98e) 

SXn-i^0 + O{A^), (98f) 

SXn = Q + 0{A^). (98g) 

As there is no scale dependent logarithmic divergence at one-loop order in three-dimensions, all the beta functions 
vanish |33[ . 

VII. DISCUSSION 

In this paper we have generalized and applied a method based on the DeWitt-Schwinger proper time expansion to 
calculate the ultraviolet divergences of the one-loop effective action associated to the RD equation. This particular 
approach involves the physical degrees of freedom plus the "ghost" fields, which are needed to account for the functional 
Jacobian that arises from a certain change of variables [H . For RDs this Jacobian is generally non- vanishing and must 
be taken into account in the computation of the characteristic functional. The importance of the effective action lies 
in the fact that it is the quantity needed to derive equations of motion, which correctly take into account fluctuations 
and interactions to a given number of loops. The effective action encodes the dynamics of the system. By contrast, 
the effective potential can only tell us about static solutions. In order to know whether the minima of the effective 
potential are relevant as solutions of the late time behavior of the system we must see how accessible these solutions 
arc. But before any of these questions can be answered, the effective action must be calculated. 

The one-loop effective action is given in terms of a functional determinant, which must be regulated and renormal- 
ized. The heat kernel technique is an established method (used in QFT) for carrying this out. In QFT the functional 
determinant appearing in the one- loop effective action is usually quadratic in both time {df) and space derivatives 



(V^). In passing to a Euclidean spacetime, the corresponding proper time equation for the kernel to be solved (29) is 
a heat equation for diffusion in n = d -I- 1 dimensions, with s playing the role of diffusion time. This is why its Green 
function (whether exactly or approximately calculated) is justifiably denoted as the heat kernel. However, for SPDEs, 
such as those considered in B , the functional determinant in the corresponding one-loop effective action involves not 
only a "mismatch" between the number of spatial and temporal derivatives, but also fourth- or even higher-order 
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spatial derivatives. We have seen an explicit example of this in the RD equation treated here, which yields second 
order in time but fourth order in space derivatives, (V^)^, for its associated operator determinant. (Even higher 
derivatives will be encountered in the one-loop effective actions based on the Sivashinski and the Swift-Hohenberg 
equations: two time derivatives but eight spatial derivatives.) While much is known about the standard heat kernel 
and its associated Seeley-DeWitt coefficients, very few (as far as the authors are aware) of these ideas have been 
applied to other types of field theories whose fluctuations may be of a non-quantum nature (i.e., noise) |29f| . The heat 
kernel technique (and its generalization) is well suited to regularize one-loop effective actions, and therefore, is very 
useful to handle theories with derivative-type interactions, as well as higher derivative "kinetic" terms. 

We have applied these techniques to compute the one-loop effective action for the RD equation. As regards its 
ultraviolet renormalizability, we found that the terms appearing at one-loop order have the same structure, i.e., involve 
the same terms present already at the level of the "classical" or zero-loop action. Strictly speaking, this claim holds 
true only if a certain bare constant is added to the original equation of motion, as we have seen. This constant, or 
"tadpole" , is needed to carry out the one-loop renormalization of the leading divergence that appears in the effective 
action. Moreover, this constant admits a simple physical interpretation and can be ascribed either to a constant flux 
rate or as the mean value of the additive noise source. In regards to the scale of application of the RD equation, 
the one-loop renormalizability indicates that although RD is a macroscopic equation, [only intended to make physical 
sense for scales greater than a certain minimum length scale Lq, defined by the underlying molecular physics (if one 
is discussing chemical reactions)], we have shown in this paper that the ultraviolet behavior of the RD equation is 
controlled, and that considered as a strictly phcnomcnological equation its short distance behavior is much better than 
one had any right to expect. The short distance structure of this effective action has been revealed via the calculation 
of the Seeley-DeWitt coefficients up to one-loop order in the noise amplitude. By means of this information, we have 
been able to establish the one-loop finiteness of RD equations driven by additive white noise for d = and d — 1 space 
dimensions and their renormalizability for d = 2 and d = 3 space dimensions. In d = 2 space dimensions there are 
logarithmic divergences which lead to running, (scale dependence) of the parameters that describe the RD equation. 
There is no wavefunction renormalization at one-loop order. These results hold for polynomial reaction kinetics of 
arbitrary degree N. (Note: The absence of wavefunction renormalization has already been demonstrated for the case 
of a (complex) RD equation describing pair reactions {N = 2), where it turns out that the one- loop beta function in 
d = 2 is exact Q.) When taken as a model for pattern development, this result becomes even more striking since this 
means we can safely use the RD equation to investigate the important short distance and small time limit of the field 
correlations that arise in pattern formation, as already remarked earlier. 

The RGE results obtained here are identical to those obtained (by different means) from an effective potential 
calculation for RD equations presented in Ref. 1^ . The effective potential is the effective action restricted to constant 
field configurations and plays an important role in uncovering patterns of symmetry breaking and in the onset of 
pattern formation. Nevertheless, the claim of the one-loop renormalizability made in ||5|] requires the investigation 
of the wavefunction renormalization which was beyond the scope of that paper. The work presented here is also 
intended to complete and complement that discussion. Moreover, as pointed out there, the combined effects of 
noise and interactions is to shift the symmetric states of the system, as well as to change the nature of the linear 
instabilities that may be induced by perturbations around these new states. A spatial pattern is, by definition, a 
spatially inhomogeneous configuration with a higher or lesser degree of symmetry, if any such symmetry is at all 
present. Thus, to investigate the onset of the spatial-temporal patterns resulting from fluctuations and interactions, 
one must go beyond the effective potential. This requires working with inhomogeneous fields and the effective action. 

The cautious reader will have noticed that most of the calculations developed here depend solely on the reaction 
potential V and its derivatives and not on the fact that V is a polynomial. In fact, this is easy to see from the 



structure of our Seeley-DeWitt coefficients, (59a)-(59b) and (60a)-(60b). Thus, the question arises if our treatment 
can be extended and applied to handle RD equations with non-polynomial reaction kinetics. The answer is in the 
affirmative provided that the potential admits a real solution to ([iq ) {e.g., V[(j)] ~ sinh[(/)]) since we need a constant 
background about which to expand the action, as indicated in (|62D. Provided such a solution exists, the rest of the 
analysis presented here carries through as is, up to the extraction of the renormalization group equations (which will 
again be non- vanishing only for d ~ 2). At this point the explicit functional form of the potential changes the nature 
of the basis set of independent field operators that leads to the RGEs. Non-polynomial reaction kinetics do indeed 
arise in many applications in chemical physics and in the modelling of biological pattern formation, typically in the 
form of rational functions (i.e., ratios of polynomials) and whenever constraints are to be imposed on the model [|l|,0. 
When a coarse-grained field theoretic approach is applied to density waves in earthquakes for example, a stochastic 
PDE for the (scalar) slip field (which measures the relative displacement of two elastic media in contact taken along 
the surface of contact) results which depends on the cosine of the slip field, cos (j>, and is driven by additive white 
noise [Q. Thus, the work presented here is broad in scope. 

Finally, these results also have the following practical application: as analytic calculations can be carried only so 
far, it is clear that numerical studies of SPDEs are crucial. Ultraviolet renormalizability corresponds to the situation 

17 



in which long distance physics is largely insensitive to the details of short distance physics. In considering numerical 
studies of the RD equation, we can therefore assert the cut-off insensitivity of the numerical solutions, at least to 
one-loop. (In numerical studies, the ultraviolet cut-off is provided by the lattice spacing.) This is of paramount 
importance since a numerical study of RD will give us the information needed to see if the system thermalizes, if it 
has steady state solutions, and most importantly, if the minima of the effective potential calculated in Ref. B are 
explored in the time evolution of the system. 
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APPENDIX A: FREE GREEN FUNCTION FOR THE RD EFFECTIVE ACTION 



In this Appendix we calculate the free Green function appearing in ( |33| ) . The only "difficult" part of the analysis 
is dealing with the fourth-order spatial bi-harmonic operator (V^)^. Using translational invariance to set x — x' —> x, 
and introducing Fourier transforms (d"fc = d'^fc dco) as follows 

/^ni. (-+00 iQ ^ ^ 

T;^^ / — Gh,4kLu\n) exp[i{k ■ X ~ Lut - Qs)], (Al) 

we easily find that 

G.ee(A^.|^) = ,.^(,,\)._,^ . (A2) 



where k^ — k ■ k. By inverting the Fourier transform (|Al| ) one obtains the following integral representation for the 
free Green function 

gfrec(a:,0|.)^/^/^/g ^,^J,^,^^^ exp[z(fc.x-u;t-».)]. (A3) 

We first perform the integral over f2 by means of a contour integral in the complex fi-plane. There is only one simple 
pole on the negative imaginary fi-axis and we close the (semi-circular) contour (centered at the origin) in the lower 
half plane. As the radius of this arc goes to infinity, only the contribution along the real fi-axis remains (s > 0). As 
a result and by the Residue Theorem we have (the contour is closed in the clockwise sense) 

J2^J h c^pW^-^-'^*)] exp(-s[c.2 -f (z.fc2)2]). (A4) 

We can go further and compute the lo integral exactly to obtain 

Gft.oc(a;,0|s) = ( - — ) exp ( -— ) / — -—T exp[-s(z^fc2)2 + ifc.fj , (A5) 



AnsJ \ is J J {2tt)'^ 

and the remaining momentum integral is manifestly convergent (for d> 0). As for the boundary condition, note that 
for s —f 0, linis^o G^*'^°°-'(xi, 00|s) — S{t,0)S'^{x,6), (in the sense of distributions) as it must, since 

and 



<5(i) ~ Inn (^— ) expHV(4.5)], (A6) 
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6'{x) 



(2^ 



(A7) 



One can also check that the boundary condition is satisfied, before integrating over w, by simply setting s = in 
(A4). Let us now work out the momentum integration. 
The angular integration is given by 



dfld-i exp lik-xj — (27r) 



.d/2 J(d-2)/2(kx) 



The exact Green function or kernel for our fi-ee "heat" operator in (p2|) is 



G,,ee(B,00|.) = ( ^) ' exp (-^) ^^ j^ 



(A8) 



k'^-^dk exp[-s(j/P)2] 



J(d-2)/2{k\x\) 



(fc|2.|)(d-2)/2 



(A9) 



This solves the differential equation (p3) and satisfies the boundary condition Gircc{x^t\Q) — S'^{x) S{t) for vanishing 
proper time s, hence this is the unique solution of (|32|). 

Important point: as remarked above, translational invariance implies that Gfrccixt,yt'\s) = Gfrcci\x — y|, (t — t')\s). 
We are treating stochastic processes on a flat d + 2-dimensional background (d-space dimensions plus real physical 
time t plus Schwinger proper time s). 
Using the Taylor series representation 



Mz) 






2n 



(-)" (^/2) 
2/ '^ n\T(i + n + l)' 



(AlO) 



and integrating (A9) term- by-term, we find that (after making use of the time and space translational invariance of 
the Green function) 



Gfrco(^i,^'t'|s) 



1 

47rs 



= Ad exp 



/ 2\-i- 
\sv ) 4 exp 


{i-m 




4s 





{t-t') 



l\2 



4s 



1 ^/ 1 

n=0 ^ 



r(f + |) /|f-x'r" 



n!r(|-|-n) V ^Vs 



/ 1 \ 



n=0 



r(f) 



n!r(f +?i) V ^> 



where 



Arf = ( ^ I — — . . ^, ,,, I ^ ) , and Cd,n 



47ry 2(27r)'*r(|) V^' 



p/ d+2rt 'i 



(All) 
(A12) 

(A13) 



Special attention should be paid to the fact that this free Green function involves a series in half-integer powers of 
Schwinger proper time, -v/s; a feature that we use in choosing our ansatz for the full Green function. 



APPENDIX B: HOMOGENEOUS FIELD CONFIGURATIONS 

For constant (homogeneous and static) fields, (f) — vq, the associated "heat kernel" can be solved for exactly. We 
simply present the final result and skip the details of a calculation that is an analog of that for Gircc{x,x'\s). The 
on-diagonal Green function Gf^^^^{x ^ x\s) is given by 

41(^.^1^) = Ad s-^--i exphs(y"[«o]l^M +rN]r[«o])]E(^2 V^'N)' ^ 

= GiUx.x\s) exp [-siV"[vo]V[vo] + V'[vo]V'[vo])] ^(^/i 2 V'Uvo])' %^, 
where we have again made use of the definition Cd/- We conclude that 



e=o 



(Bl) 
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+ 00 +00 „ 

^[6|] //2 = exp [-S {V"[v,]V[vo] + V'[vo]V'[vo])] ^(Vi 2 V'MY ^. (B2) 

£=0 1=0 

We now match the first fractional powers in s and obtain the Seeley-DeWitt coefficients for a constant field configu- 
ration Vq 

[bo] = 1, (B3) 

[bi] = 2Cd,iV'[vo], (B4) 

[bi] = -V"[vo]V[vo]. (B5) 

The coefficients for G^-^^^^. can be immediately obtained from the previous Seeley-DeWitt coefficients by setting 
V"[vo] = 0. That is 

[ao] = 1, (B6) 

[ai] = 2 Cd,i V'[vol (B7) 

K]=(^-l) V'[vo]V'[vo]. (B8) 

Finally, for the net Seeley-DeWitt coefficients, we can write 

[Co] = 0, (B9) 

[ci] = 0, (BIO) 

[ci] = -V"[vo] V[vo]. (Bll) 

These results are consistent with the net integrated Seeley-DeWitt coefficients an d w it h th e physical and ghost 
integrated Seeley-DeWitt coefficients presented in the body of the paper [equations ( |59a| )-( |60q )]. 

APPENDIX C: THE FEYNMAN-HELLMAN FORMULA 

We write the Feynman-Hellman formula in the form [pOJ 

— exp(^ + eB)= / de exp[^(A + eB)] B exp[(l - i){A + eB)]. (CI) 

de Jo 

This equation is central to the computation of the Seeley-DeWitt coefficients presented in the paper. For instance, if 
we take the trace and then use the cyclic property, we can write 

— Tr[exp(A + eB)] = Tr {B cxpUA + eB)]} . (C2) 

de 

We differentiate the previous equation to obtain 



de2 



Tr[exp(A + eB)] ^ f dt Tt{B exp[i{A + eB)] B exp[(l - £){A + eB)]} . 
Jo 



(C3) 



We can conclude that 

^2 (-1 



Tr[exp(A + eB)] = Tr[exp(A)] + eTr[B exp(yl)] + — / di Tr {B exp(M) B exp[(l - e)A])} + 0{e^). (C4) 

2 Jo 

This perturbative expansion is the basis for extracting the first two integrated Seeley-DeWitt coefficients. 
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